Long range order in the classical kagome antiferromagnet: effective Hamiltonian approach 
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Following Huse and Rutenberg [Phys. Rev. B 45, 7536 (1992)], I argue the classical Heisenberg antiferro- 
magnet on the kagome lattice has long-range spin order of the i/3 x a/3 type in the limit of zero temperature. 
I start from the effective quartic Hamiltonian for the soft (out of plane) spin fluctuation modes, and treat as 
a perturbation those terms which depend on the discrete coplanar state. Soft mode expectations become the 
coefficients of a discrete effective Hamiltonian, which (after a coarse graining) has the sign favoring a locking 
transition in the interface representation of the discrete model. 

PACS numbers: 75.10.Hk, 75.50.Ee, 75.40.Cx 



I. INTRODUCTION 

Consider the nearest-neighbor antiferromagnet with clas- 
sical spins of n = 3 components on the kagome lattice of 
corner-sharing triangles, 



H 



(ij) 



This is the prototypical highly frustrated system, meaning its 
ground state manifold has macroscopically many degrees of 
freedom, and any spin order or freezing sets in at tempera- 
tures T < J-. It is well established that as T — > 0, the free 
energy of spin-mode fluctuations causes ordering into a copla- 
nar state, a particular kind of classical ground state in which 
all spins lie in the same plane of spin space pointing in just 
three directions (120° apart) 2 . These directions - which can 
be written as colors Cj = A, B, or C taken by spins in a 3-state 
Potts model - constitute a 3-coloring (the ground state con- 
straint implies every triangle has one of each color). The num- 
ber of such colorings is exponential in the system size. The 
same is true for three-dimensional lattices of corner-sharing 
triangles such as the (half) garnet lattice 3 or equivalently hy- 
perkagome lattice 4 , and others^. 

Can the coloring achieve a long-range order? All simu- 
lations 2 ' 6 ' 8 ' 9 indicate the Potts spins are disordered (or alge- 
braically correlated) as in the unweighted coloring (see be- 
low). However, following Huse and Rutenberg 10 , I propose 
this coloring develops long-range order in the T — > limit, as 
a consequence of the unequal weighting of the discrete states, 
when one takes into account the free energy of fluctuations 
about each state. Of course, for d = 2 at T > 0, the orien- 
tation of the spin plane must fluctuate slowly in space; never- 
theless the colors/Potts directions may be unambiguously de- 
fined throughout the system. But my goal is only the T — > 
limiting ensemble, well defined (and nontrivial) since the ob- 
tained effective Hamiltonian [e.g. ( [ToT ll scales as T; whereas 
the spin-plane correlation length diverges exponentially-^ as 
T -> 0. 

The calculation entails a series of mappings and effective 
Hamiltonians. First I shall review how, starting from the usual 
spin-deviation expansion, one integrates out most of the fluc- 
tuations leaving a quartic effective Hamiltonian Q for the 
dominant fluctuations 11 . The largest terms of Q are indepen- 
dent of the discrete Potts configuration, so treating the rest as 



a perturbation yields an effective Hamiltonian $ for the Potts 
spins, purely entropic in that $ oc T. Its coefficients may be 
inferred from simulations, or approximated analytically (tak- 
ing advantage of a "divergence constraint" on the dominant 
fluctuations). The Potts spins map in turn to a "height model", 
whence it becomes clear that $ causes locking into an ordered 



n\ stated*. The expected long-range order is too tenuous to see di- 
rectly in simulations, but might be estimated analytically from 
the height model. 



II. EFFECTIVE HAMILTONIAN DERIVATION 

The object is to obtain an effective Hamiltonian for any of 
the discrete coplanar ground states, which absorbs the free en- 
ergy of the low-temperature (anharmonic) fluctuations about 
that state. The first step is the "spin-wave" expansion in de- 
viations from a given coplanar state. We parametrize the out- 
of -plane deviation as cr,, and the other deviation component 
as 9i, the spin's in-plane rotation about the plane normal axis. 
Then the spin- wave expansion (I set J = 1) is 

n sw = n {2) + h {3) + n {i) + . . . , 

n(2) - E[i^-^) 2 +H + 5> 2; (2a) 
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In Eq. dZbb . a indexes the center of each triangle, and 
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(3) 



From here on, I use "am" (m = 1, 2, 3) to denote the site on 
triangle a in sublattice m, as an alias for the site index "i" ; the 
index m is taken modulo 3 (in expressions like "m + 1") and 
runs counterclockwise around the triangles whose centers are 
even sites on the honeycomb lattice of triangle centers. Fol- 
lowing Ref. 11,1 retained only "dominant" anharmonic terms 
H^dom an d T^ 4 "* dom, being the parts of (O and d2cl ) contain- 
ing the highest powers of a (this will be justified shortly). 



2 



The r] a prefactor in H ^ is the only dependence in Eq. [2] 
on the coloring state; this "chirality" rj a , is defined by rj a = 
+1(— 1) when the Potts labels are ordered as ABC (CBA) as 
one walks counterclockwise about triangle a. It is convenient 
to label coplanar states by the configuration {rj a } 12 . Then a 
discrete Hamiltonian $ can be defined for colorings, a func- 
tion of the {ci] implicitly through the r) a 's in (0: 

e -*({ Cl }/T = z{{ci}) = I JJ ( ^. daJe-H-ZT (4 ) 

J basin ^ 

As T — > 0, the ensemble weight concentrates closer and 
closer to the coplanar state 213 ; the integral in (|4]i is implic- 
itly limited to the "basin" in configuration space centered on 
one coplanar state, and Z({ci}) is the portion of the total par- 
tition function assigned to the corresponding coloring. Since 
TL^> is independent of which coplanar state we are in, $ is 
independent of {a} at harmonic order. 

Before we go on to anharmonic order, let's note the a part 
of W< J ) can be written = \ £ Q (£ m o am f. So there 

is a well-known whole branch of out-of-plane (a) modes, 
called "soft modes", having zero cost at harmonic order; the 
soft mode subspace is defined by the constraint 

3 

= (Soft) (5) 

m— 1 

being satisfied on every triangle a. (Two more out-of-plane 
branches, as well as all 9 branches, are called "ordinary" 
modes.) Being limited only by higher order terms, soft modes 
have large mean-square fluctuations, of O(VT), compared to 
0(T) for ordinary modes 211 ; this explains why factors con- 
taining soft modes were "dominant" in Eq. ©. The <7j's in 
"dominant" terms are limited to the "soft" subspace satisfying 
©. 

The next step is to do the Gaussian integral over all 9{ 
modes 2 , as worked out in Ref.Qjl obtaining a quartic effective 
Hamiltonian Q for only soft modes: 

Q = H {4) dom -J2v a VpQ' al 3 (6) 

with 14 

Q'a,0 = GamfinV 1 am o'}l n ( 7 ) 

m.n—1 

The Green's function of the 8 modes was defined by 

T G am ^n = ((#a,m+l _ »a,m-l) (&0,n+l ~~ "/3,n-l))e (8) 

where "(...)„" means taken in the (Gaussian) ensemble of 
H (2) e (= the 9 part of ri {2) ). As G i} decays with dis- 
tance, the largest terms are state-independent: Q' aa — 

(3/16)[GoEm=l (7 am + 2G 1 Em<n °"am°"an]- where G 

and G\ are the on-site and first-neighbor G L j. Trivially 
2Gi = —Go, and Gq = 1 (due to equipartition, which 
implies (H^e) = 3T/4 per triangle). Also, given ©, 



E ra <„ ^lEm cr«m in Q' aa > and similarly in © 

7~t^dom — * (1/16) J2i a t- Finally we can regroup © as 

Q=Qo-J2 VaVpQ^p, Qo = B J2 °i ( 9 ) 

with B Q = 13/16 from both H (4) do m and Q' aa terms. 

Now I turn to the perturbation expansion: the key step in 
our whole derivation is to expand (01 treating the {r) a } as if 
they were small quantities. (In fact \rj a \ — 1, so a pertur- 
bative treatment might appear questionable, but quantitatively 
Qo has a much larger coefficient than the terms in Q', owing 
to the decay of Gy with separation.) The resulting (and final) 
effective Hamiltonian is, to lowest order, 

1 \ - 

$ = JotpVaVP (10) 

with 

m,n— 1 

where the expectation is taken in the ensemble of Q . No- 
tice that since Q is homogeneous in {ci}, it follows that the 
partial partition function Z({q}) in (0]i - and consequently 
$ jT - is temperature independent as T — > 0, apart from a 
configuration-independent powers of T. 

A corollary of my assumption that Q! a a is "small" is that 
expectations (•••) sw of polynomials in {at}, measured under 
the full spin-wave Hamiltonian 7i sw , should be practically in- 
dependent of the coloring configuration {rj a }^. That can be 
checked in Monte Carlo or molecular dynamics simulations 17 
of the Heisenberg model. The needed correlations can be 
measured even if the system is confined to the "basin" of one 
coplanar state: there is no need to equilibrate the relative oc- 
cupation of different basins. Those same simulations would 
numerically evaluate the quartic expectations 

III. SELF-CONSISTENT APPROXIMATION FOR 
COUPLINGS AND ASYMPTOTIC BEHAVIOR 

An alternative to simulation is to analytically evaluate the 
quartic expectations in ( fTTb . using a self-consistent decou- 
pling. That is, (O is replaced by 

F val ^^Bj2^l (12) 

i 

defining a Gaussian variational approximation to the soft 
mode ensemble; here 

B ee 6B (<if) vaT , (13) 

with "(•••) var " taken in the ensemble of ( fT~2T >. Now let Ty (also 
written Y a m,pn) be the Green's function for a.- L modes: 

(o i cr j ) vax =TT ij /B (14) 
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(this definition makes independent of B and T) and let 
T u = T = 1/3. Combining (Qj) and (O, I get the self- 
consistency condition B = (QBqTqT) 1 ' 2 = (13T/8) 1 / 2 . 
Next, the expectations in (fTTb are evaluated in the variational 
approximation, decoupling by Wick's theorem as 



(°i )var(°j)var 



-2((TjCT 



(§)> 



J I var 



■ 2 ^J 

(15) 

Substituting ( fT5l ) into dTTb gives my central result for the ef- 
fective Hamiltonian, 



Ja/3 

T 



3 



m,n— 1 



(16) 



[Tg from (TTSb always cancels in the m, n sum.] 

Eq. (O gives Ji/T « -1.88 x lO" 3 and J~ 2 /T » 
—4.3 x 10~ 4 . Assuming these two dominate, the state with 
the lowest $ value is the "V3 x v3" pattern, the "antifer- 
romagnetic" arrangement of chiralities rj a , but other coloring 
configurations are only slightly less likely. The long range or- 
der suggested by Ref. 10 is a subtle crossover of correlation 
functions at large (but not diverging) scales, best expressed in 
terms of a "height model", as will be developed in Sec.lIVI 

Before that, in order to check that more distant couplings 
J a j3 can be neglected, I will work out how they scaling at large 
R. We need both kinds of Greens function in (fl6t . tackling the 
Qi fluctuations first. In Eqs. (|2ab and (0, (# m +i — # m -i) ~ 
—a e a e^ • V 8, where a is the nearest neighbor distance, and 
e a = +1(— 1) when a labels an even (odd) triangle. The unit 
vector e m = (cos^mj smip m ), is defined to point from the 
center of any even triangle to its m corner, and e„ = i x e m . 
At long wavelengths, 



i 



d 2 r\Vd(r) 



(17) 



where pg = \[Zj2. Asymptotically the Greens function of 
(fTTT i is pseudo-dipolar: 



G 



am,{3n 



2irp 6 R 2 



e a e/3 cos(ip m + ip n - 2i/j r ). 



(18) 



Here (R, V>_r) are the polar coordinates of the vector between 
triangle centers a and (3. 

The Gi fluctuations are handled similarly. The soft-mode 
constraint (0 is implemented by writing er s ; as a discrete 
gradient, <Tj = <j) v — (p^, analogous to the "height" model 
constructions 1 - Here {</> M } is defined on the hexagon cen- 
ters, and [i — > v is oriented counter-clockwise around even 
kagome triangles. The discrete gradient defining <x; can be 
converted into a continuous one, a am « 2a ■ V0. Then 
the long-wavelength limit of (fT2b is looks like dTTb . with 
Pe — * P4> = 2^/3 B. That implies that for large separations 
R, ^am,/3n looks like Eq. ( fT8l with pg — > p$/4:. Inserting 
both Green's function behaviors into ( fT6] i, I get the asymp- 
totic behavior of the couplings: 



a/ 3 



(19) 



for large R with A = &^/1??tt 3 w 2.0 x 10" 3 . Eq. © 
shows the interaction decays rapidly with distance and oscil- 
lates as a function of angle. 



IV. HEIGHT MODEL AND LONG RANGE ORDER 

The discrete ensemble in which all 3-colorings {ci} are 
equally likely is known to have power-law correlations, which 
may be understood via a mapping of the Potts microstates to a 
two component "height" variable h(r) 10,19 . At coarse-grained 
scales, the ensemble weight of {h(r)} is described by a free 
energy 



F h = I d 2 v±K\Vh\ 2 , 



(20) 



handled by standard Coulomb-gas techniques^- 

Ref. LL0| pointed out the equal-weighted coloring has a 
height stiffness K = K c exactly, where K c is the critical 
value for the roughening transition. Any increase in K must 
cause h(r) to lock to a uniform mean value. 1 0,20 ' 21 . That cor- 
responds to long-range order of the colors (= Potts spins), into 
the pattern of with the flattest h(r), namely the "vo x V3" 
state. Since (as shown above) $ favors that flat state, the col- 
oring ensemble with the $ weighting is coarse-grained to a 
height ensemble with a slightly larger K, and therefore we 
get long range order, as claimed^ 

The couplings J a p as approximated analytically, or ob- 
tained from a simulation, may be used as a Hamiltonian in 
discrete simulations of the coloring model. These are far faster 
than simulations of the Heisenberg spins, but I still doubt such 
simulations will see long-range order directly, in the accessi- 
ble system sizes. But the height stiffness K can be accurately 
measured (using Fourier transforms 22 .) With that, by iteration 
of renormalization-group equations 21 , it should be possible to 
semi-analytically estimate the length scale £ at which the color 
correlations cross over from power-law decay to long-range 
order, and the size of the order parameter. 

What happens to this whole story in d — 3, for the Heisen- 
berg antiferromagnet on triangle-sharing lattice s 3 ' 4 ! 5 ? A minor 
difference is that in d — 3 the spin plane orientation has true 
long-range coplanar order at some T > 0, as do the three 
spin directions within the plane 8 . The derivation and result 
for the effective Hamiltonian ( [ToT l extend to d = 3; There is 
also the unimportant difference is that, in deriving the asymp- 
totic behavior of J a p, a "Coulomb phase"^2i rather than a 
"height function" viewpoint must be used for coarse-graining 
a, but Tij still has a pseudodipolar for m 23 ' 24 , and the final 
asymptotic form is analogous to ( fT9l (J a p oc 1/R 9 with an 
oscillating angular dependence). 

The crucial difference in d = 3 is that the discrete (Potts) 
variables also have a "Coulomb phase" in place of the "height 
representation" used by Ref.[10|. There exists a coarse-grained 
"flux field" analogous to Vh, but the analog of h itself is a 
vector potential and is not uniquely defined. The Hamiltonian 
$, I conjecture, tends to favor states with zero coarse-grained 
flux, which means it tends to increase the flux stiffness K of 
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the three-dimensional model. But, in contrast to two dimen- 
sions, in the absence of $ the system is not sitting at a critical 
K; therefore, the tiny increase in K due to the transverse spin 
fluctuations cannot drive us into a new phase. Thus, no long- 
range order of the colorings is expected in d = 3, merely the 
the pseudodipolar correlations inherent to the Coulomb phase. 

V. CONCLUSION 

A path has been shown to the elusive long-range order 
of the classical kagome antiferromagnet, through a string of 
mappings or elimination of degrees of freedom: ground states 
to colorings to chiralities to discrete h M height representation 
and finally its coarse-grained continuum version. Other maps 
go from all spin deviations, to soft modes Cj, to their height 
field (f)^ or 4>(y). The boldest approximations were (i) the 
perturbation expansion (ITTb of the effective quartic Hamil- 
tonian ©; this had no controllable small parameter, but it 
was argued the terms were numerically small (ii) the varia- 
tional/decoupling handling of the quartic ensemble {&i}. In 



place of the approximations used here, the more elaborate but 
more controlled large- TV approach 24 ' 26 (where N is the num- 
ber of classical spin components) looks promising as a formal 
way to vindicate both approximations. 

The philosophy followed here 2 ^ is to obtain an effective 
Hamiltonian defined for arbitrary spin arrangements, not just 
specially symmetric ones (even if that necessitates cruder ap- 
proximations). I have previously used the trick of turning the 
spin configuration into a set of coefficients or matrix entries 
and then expanding in them for several system s 15,27,28 . In par- 
ticular, a related expansion in Ti 1 - 3 ^ to obtain a Hamiltonian of 
form ( [Tol l was carried out for the large-S* quantum Heisenberg 
antiferromagnet in Ref . l28i 
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